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ACOUSTIC  PROPAGATION  IN  A  WEDGE-SHAPED  OCEAN  WITH 
PERFECTLY  REFLECTING  BOUNDARIES 


INTRODUCTION 

The  ocean  overlying  the  continental  slope  or  a  sloping  beach  is  most  simply  represented  as  a 
wedge-shaped  domain  with  perfectly  reflecting  (pressure-release)  boundaries.  This  model  is,  of  course, 
inadequate  in  that  it  cannot  account  for  certain  phenomena,  such  as  acoustic  penetration  of  the  bottom, 
which  are  encountered  in  the  real  ocean.  On  the  other  hand,  the  wave  equation  for  the  perfect  wedge 
is  separable,  and  the  analysis  of  the  acoustic  field  due  to  a  point  source  within  the  domain  is  tractable. 
By  way  of  contrast,  the  wave  equation  for  the  penetrable  wedge  is  not  separable,  a  fact  which  inevitably 
complicates  the  analysis  of  the  acoustic  field  in  this  case.  Here  we  shall  confine  our  attention  to  the 
simpler  situation  where  the  boundaries  are  perfect  reflectors.  Certain  physical  attributes  of  the  perfect 
wedge  are  discussed,  most  notably  the  fact  that  the  radiation  field  is  constrained  to  form  modal  beams 
as  a  result  of  the  multiple  acoustic  interactions  with  the  boundaries.  This  conclusion,  which  derives 
from  the  solution  of  the  wave  equation,  can  be  most  satisfactorily  interpreted  by  calling  on  the  concept 
of  ray /  mode  duality.  The  question  of  the  connection  between  rays  and  modes  in  the  wedge  is  pursued 
here  at  some  length. 

In  general,  the  radiation  field  in  a  wedge-shaped  domain  formed  by  two  perfectly  reflecting  bound¬ 
aries  consists  of  two  components:  the  modal  component  is  a  consequence  of  the  constraints  imposed  by 
the  boundary  planes  on  the  field,  and  the  diffracted  component  arises  from  scattering  at  the  apex  of  the 
wedge.  It  is  well  known  [1]  that  when  the  wedge  angle  is  a  submultiple  of  it,  the  diffracted  component 
is  identically  zero;  and  in  certain  special  cases,  even  when  the  wedge  angle  is  not  a  submultiple  of  n, 
the  diffracted  component  is  again  absent  [2],  What  is  more  significant,  however,  is  that  in  the  context 
of  ocean  wedges,  where  the  wedge  angle  is  usually  very  small  (typically  a  few  degrees)  and  the  ranges 
of  the  source  and  the  receiver  from  the  apex  are  generally  many  wavelengths  long,  the  diffracted  com¬ 
ponent,  if  not  absent,  is  negligibly  small.  This  being  so,  it  is  ignored  in  the  following  discussion,  which 
is  devoted  to  the  modal  component  of  the  field. 

An  approximate  solution  for  the  modal  part  of  the  field  in  a  perfect  wedge  has  been  derived  by 
Bradley  and  Hudimac  [3]  and  examined  by  Graves  et  al.  [4],  Unfortunately,  its  range  of  validity  is 
extremely  limited.  In  particular,  it  is  unsatisfactory  throughout  a  large  proportion  of  the  modal  beams 
alluded  to  above,  in  just  the  regions  which  are  of  interest  in  many  applications.  The  failure  of  Bradley 
and  Hudimac’s  solution  in  these  regions  could,  perhaps,  have  been  anticipated,  since  it  does  not  show 
the  required  form  around  the  source  point.  One  of  the  main  results  in  the  present  report  is  a  new  solu¬ 
tion  for  the  field  in  the  wedge,  which  has  a  much  more  extensive  range  of  validity  than  that  of  Bradley 
and  Hudimac,  and  which  behaves  correctly  in  the  immediate  vicinity  of  the  source. 

THE  FIELD  DUE  TO  A  POINT  SOURCE  IN  THE  WEDGE 

The  geometry  for  the  wedge  problem  is  shown  in  Fig.  1,  where  S  and  R  represent  a  point  source 
and  a  point  receiver,  respectively.  Since  we  are  interested  in  the  three-dimensional  field  in  the  wedge, 
S  and  R  are  not  necessarily  in  the  same  vertical  plane  running  perpendicular  to  the  apex  of  the  wedge. 
The  angle  of  the  wedge  is  0O.  A  cylindrical  coordinate  system  is  used  for  the  problem,  with  the  z  axis 
running  down  the  apex  of  the  wedge  and  the  angular  coordinate  9  measured  from  the  surface.  Thus 
the  boundary  planes  are  at  6  =  0  and  9  =  60.  The  coordinates  of  the  receiver  are  ( r,6,z )  and  those  of 
the  source,  which  for  convenience  is  located  in  the  z  =  0  plane,  are  the  primed  quantities  (r',9',0). 
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APEX  IS  THE  SHORE  LINE 


S  —  source 
R  —  receiver 


Fig.  I  —  Cylindrical  coordinate  system  for  the  wedge  problem 


The  equation  that  must  be  solved  for  the  field  in  the  wedge  is  the  inhomogeneous  Helmholtz 
equation: 

V2<f>  +  k20  =  -  Q—  ~r'}  8 (9  -0')8(z),  (1) 


where  Q  is  the  source  strength,  8()  is  the  Dirac  delta  function,  and  k  is  the  wavenumber  (i.e., 
k  =  (o/c,  where  «  is  the  angular  frequency  and  c  is  the  speed  of  sound  in  the  medium,  which  is 
assumed  to  be  independent  of  position).  The  function  0  in  Eq.  (1)  is  the  velocity  potential  excluding 
the  time-dependent  factor  in  the  case  of  an  harmonic  source,  or  it  is  the  Fourier  transform  (with 
respect  to  time)  of  the  velocity  potential  for  an  impulsive  source.  On  expressing  the  Laplacian  operator 
in  cylindrical  coordinates,  Eq.  (1)  becomes 


I  JL 

r  Br 


rB$ 


Br 


+  ^^  +  ^  +  k2<l>  =  ~Q  ~  ~  8(6  ~  e')8iz)- 
r  B9 1  Bz L  r 


(2) 


Equation  (2)  is  separable  for  the  perfect  wedge,  and  it  can  be  solved  for  <f>  in  several  ways.  The 
method  used  here  is  to  apply  a  sequence  of  integral  transforms  to  both  sides,  and  subsequently  to  apply 
the  inverse  transforms  to  arrive  at  the  solution.  Assuming  pressure-release  boundaries,  the  first  of 
these  transforms  is  the  finite  Fourier  sine  transform: 

0s  =  J '0°  v  OdQ,  (3a) 

whose  inverse  is 


4>(9)  =  —  V  (v)  sin  v  9. 
Oo 


(3  b) 


In  these  expressions,  v  =  m7r/0o,  where  m  is  a  positive  integer,  and  the  symbol  E„  means  a  sum  over 
all  possible  values  of  m.  It  is  clear  from  the  structure  of  Eq.  (3b)  that  the  final  solution  for  the  field 
will  be  a  sum  of  normal  modes,  and  that  m  is  the  mode  number.  Note  also  that  this  formulation 
ensures  that  the  pressure-release  boundary  conditions  are  satisfied.  On  taking  the  sine  transform  of  Eq. 
(2),  we  obtain 


B_ 

Br 


Br 


vl  ,  ,  9_0s 
r2  Bz2 


+  k2(f>s  =  —  Q  — —S(z)  sin  v  9'. 


(4) 


The  next  transform  to  apply  is  the  generalized  Hankel  transform  of  order  v  =  mn/90,  which  we 
shall  assume  to  be  an  integer.  This  is  tantamount  to  saying  that  90  is  a  submultiple  of  tt,  so  that  the 
field  due  to  diffraction  at  the  apex  of  the  wedge  is  zero.  The  Hankel  transform  is  defined  as 
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Ji  OO 

Q  r  <f>s(r)  J„(pr)dr, 

and  the  inverse  transform  is 

J.  oo 

0  P  4>v,s  (p) h (pr)dp. 

An  important  property  of  the  Hankel  transform  is  that  it  acts  on  the  quantity  A „<f>s  as  follows  [5]: 


where  A„  is  the  differential  operator 


=  -p2hms), 


A  = 


i  d 

d 

_ 

r  dr 

"  dr 

r 

Thus,  when  Eq.  (4)  is  Hankel  transformed  we  find  that 

d2<t>v, 


dz2 


+  ( k 2  —  p2)  4>v  s  =  —  Q  sin  v  9'  Jv(pr')d(z). 


(5a) 

(5b) 

(6a) 

(6b) 

(7) 


The  final  transformation  to  apply  is  the  Laplace  transform,  defined  as 

Ji  oo 

<£„,j (z)  exp  ( -qz)dz , 


(8) 


where  the  unilateral  form  has  been  chosen  because  the  field  must  be  symmetrical  in  z.  On  applying  this 
integral  transformation  to  Eq.  (7),  we  find  that 


$/,  v,s 


q<t>v  s  (0)  — ^  sin  v  0 '(/„(/?/•') 


(9) 


( q 2  +  k2  —  p2) 

where  0)  is  the  value  of  <f>v  s  at  z  =  0.  Now,  the  inverse  Laplace  transform  of  Eq.  (9)  is  easily 
shown  to  be 


tfv.s  =  tfL.stO)  cos  7)|z|  -  sin  v  9'  Jv{pr')  —  T?'Z|  , 

2  T) 


where 


V 


=  Vk2  —  »2. 


(10) 

(ID 


From  the  radiation  condition,  that  is,  the  requirement  that  the  field  should  go  to  zero  when  \z  \  — ►  °o 
for  p>k ,  we  find  that 


and  hence 


0„,s( 0)  =  j  sin  v  9'  ~--Pr  ^ , 
2  7) 


<t>v,s  =  j  0/2  sin  v  0'  -  exp  (j  tjIzI). 


(12) 


(13) 


The  next  inversion  integral  to  apply  is  the  inverse  Hankel  transform  defined  in  Eq.  (5b).  Using 
the  expression  in  Eq.  (13),  this  gives 


O  c  ‘ 

4>s  =  J  j  sin  v  6'  Jo 


p  exp  {j7]\z\)Jv{pr)Jv{pr') 


dp. 


(14) 


Finally,  we  obtain  the  solution  for  the  field  in  the  wedge  by  taking  the  inverse  sine  transform  of  <J>S 
defined  in  Eq.  (3b): 


Y  I„(r,r' ,z)  sin  v  9  sin  v  9'. 
tin 


(15) 
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This  is  a  sum  of  normal  modes  whose  coefficients  are  given  by  the  integral 

I v(r,r',z)  =  j  JQ  p  -Xp  Jv(pr)J„(pr')  dp.  (16) 

Note  that  the  solution  as  it  stands  in  Eq.  (15)  is  exact,  since  no  approximations  have  yet  been  intro¬ 
duced.  It  is  only  in  evaluating  Iv  that  some  accuracy  will  be  sacrificed. 


THE  NORMAL  MODE  INTEGRAL 

To  evaluate  the  integral  in  Eq.  (16),  it  is  convenient  to  convert  it  to  another  form  using  the 
Bessel  function  identity  [6] 

J„(pr)J„(pr')  =■*  —  f  J0(p  Vr2  +  r'2  —  2 rr'  cos  cr)  cos  v  a  da,  (17) 

77 

which  is  valid  for  v  an  integer.  When  this  expression  is  substituted  into  Eq.  (16),  we  obtain  the  double 
integral 

/„  =  f  cos  v  a  J  p  6XP  j  (p  Vr2  +  r'2  —  2rr'  cos  cr)  dpda.  (18) 

77  *'0  *,0  T) 

Now,  the  inner  integral  here  is  the  field  due  to  a  point  source  in  an  infinite  medium  [7],  and  hence  /„ 
can  be  written  as 


where 


l  r77 

nRo  *'0 


cos  v  a 


exp[— jkRo(l  —  2 a  cos  cr) 1/2] 
(1  —  2 a  cos  c r)1/2 


da, 


(19) 


R0  =  (r2  +  r'2  +  z2)1/2,  a-  —  ^  1/2.  (20) 

From  the  definitions  in  Eq.  (20),  it  is  clear  that  a  takes  its  maximum  value,  equal  to  1/2,  when  r  =  r* 
and  z  =  0.  As  the  source  and  receiver  separate  in  the  r  —  z  plane,  a  falls  in  value,  approaching  zero  in 
the  limit. 


The  form  of  the  integral  in  Eq.  (19)  is  more  amenable  to  evaluation  by  approximation  techniques 
than  is  the  form  in  Eq.  (16).  An  evaluation  procedure  is  described  in  Appendix  A  which  leads  to  the 
result 


J  1  .VI T 

h  -  j  exp  j— 


exp  1-jfcR^l  +  b\)}H^  {kR\b\) 


Ri 


+ 


exp  {-jkR2 (1  -  b2)}H ®  ( kR2b2 ) 


R 7 


(21a) 


which  is  valid  provided  the  condition 


kR0a  »  1  (21b) 

is  satisfied.  In  the  expression  for  /„  in  Eq.  (21a),  ( )  and  //„(2)  ( )  are  Hankel  functions  of  order  v 

of  the  first  and  second  kind,  respectively,  and  the  parameters  /?,  and  bjt  i  =  1,  2  are  defined  as  follows: 


R\  =  *0(1  -  2 a)x'\  R2  =  R0(l  +  2a)1/2,  bx  =  h  =  (1  +  2a)  '  (22) 

Notice  that,  apart  from  v,  the  formulation  of  /„  given  above  depends  on  only  two  independent  vari¬ 
ables,  namely  kRo  and  a.  Equations  (15)  and  (21)  specify  the  modal  field  in  the  wedge. 
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In  the  context  of  ocean  wedges,  where  the  wedge  angle  is  usually  substantially  less  than  10°,  the 
inequality  in  Eq.  (21b)  does  not  represent  a  serious  practical  restriction  on  the  validity  of  the  solution 
for  /„.  As  we  shall  demonstrate,  each  mode  in  the  wedge  takes  the  form  of  a  beam  of  radiation  which 
diverges  as  the  range  from  the  apex  increases.  Within,  and  to  some  extent  beyond,  a  modal  beam,  the 
condition  in  Eq.  (21b)  is  always  satisfied.  It  is  only  deep  within  those  regions  which  are  not  illuminated 
by  the  source,  that  is  to  say,  the  acoustic  shadow  zones,  that,  Eq.  (21b)  may  not  hold,  but  this  is  not  a 
significant  limitation  because  the  field  in  these  regions  is  essentially  zero  anyway.  We  conclude,  there¬ 
fore,  that  Eq.  (21a)  is  a  valid  description  of  the  modes  for  most  scenarios  appropriate  to  the  ocean 
wedge. 


It  is  worth  noting  in  particular  that  Eq.  (21a)  is  a  valid  representation  of  the  field  in  the  immediate 
vicinity  of  the  source  point  (assuming  that  the  source  is  many  wavelengths  from  the  apex  so  that  equa¬ 
tion  (21b)  is  satisfied).  Under  these  circumstances  it  can  be  shown  that  the  mode  sum  for  the  field  in 
the  wedge  reduces  to  the  free-field  expression  for  a  point  source  (see  Appendix  B).  This  is  exactly  the 
form  of  limiting  behavior  that  is  required  of  the  solution. 


It  is  interesting  to  examine  Eq.  (21a)  in  another  limit,  in  this  case  as  a  — 1 ►  0,  corresponding  to  a 
large  separation  between  the  source  and  the  receiver  in  the  r  —  z  plane.  In  this  limit,  we  have 
R\  =  R2  =  Rq  and  b\  =  b2  =  a,  which,  when  substituted  into  Eq.  (21a)  give 


lim  /„  = 

o-0 


1 

Ro 


exp 


exp  -(JkR0)Jp(kR0a), 


(23) 


where  the  sum  of  the  two  Hankel  functions  has  been  expressed  as  a  Bessel  function  of  the  first  kind. 
This  expression  is  the  solution  for  the  modal  coefficients  that  was  obtained  by  Bradley  and  Hudimac  [3], 
Naturally,  its  range  of  validity  is  considerably  less  than  that  of  the  solution  in  Eq.  (21a),  since  the  addi¬ 
tional  constraint  a  =  0  has  now  to  be  satisfied.  This  condition  means  that  Eq.  (23)  is  not  a  valid 
representation  of  the  field  in  the  vicinity  of  the  source,  as  exemplified  by  the  fact  that  it  does  not  show 
the  required  singularity  at  the  source  point. 


The  implications  of  the  condition  a  =  0  are  well  illustrated  by  an  example.  If  we  assume  that  a 
<0.1,  and  that  the  source  and  the  receiver  are  in  the  same  vertical  plane  perpendicular  to  the  apex  so 
that  z  =  0,  then  a  simple  calculation  shows  that  Eq.  (23)  is  valid  only  when  r  <  r'/lO  or  when 
r  ^  10r'.  Taking  the  range  of  the  source  from  the  apex  as  1  km,  say,  this  means  that  Eq.  (23)  is 
invalid  for  all  receiver  ranges  (measured  from  the  apex)  between  100  m  and  10  km.  It  is  unfortunate, 
in  view  of  the  relative  simplicity  of  Eq.  (23),  that  this  is  precisely  the  range  that  is  likely  to  be  of 
interest  in  practical  situations. 

An  impression  of  the  relative  accuracies  of  Eqs.  (21a)  and  (23)  is  conveyed  by  Table  1,  which 
shows  the  real  part,  the  imaginary  part,  and  the  modulus  of  the  quantity  ttR0Ip  for  various  values  of  a 
and  kR0,  but  chosen  so  that  the  product  kR0a  =  100.  For  comparison,  a  numerically  computed 
evaluation  of  the  integral  in  Eq.  (19)  is  included  in  the  table.  The  agreement  between  Eq.  (21a)  and 
the  computed  values  is  apparent.  Equation  (23),  on  the  other  hand,  performs  poorly,  especially  with 
regard  to  phase.  This  is  true  even  when  a  takes  the  relatively  low  value  of  0.1  used  in  the  example 
cited  above.  An  improvement  occurs  only  when  a  is  substantially  less  than  0.1,  confirming  that  Eq. 
(23)  has  an  extremely  limited  range  of  validity. 


THE  EIGENFUNCTIONS  FOR  THE  WEDGE 

The  trigonometric  functions  in  Eq.  (15)  contain  the  entire  angular  dependence  of  the  modes  in 
the  perfect  wedge.  These  oscillatory  functions  are  analogous  to  the  eigenfunctions  obtained  for  a 
shallow-water  channel  with  pressure-release  boundaries.  In  the  shallow  water  case,  the  variable  is  the 
depth  of  the  sensor  or  the  source,  normalized  to  the  channel  depth,  whereas  in  the  wedge  it  is  the 
angular  depth  normalized  to  the  wedge  angle. 


5 


M.  J.  BUCKINGHAM 


Table  1  —  The  Modal  Coefficients  for  v  =  18  (corresponding 
to  the  first  mode  in  a  10°  wedge) 


a 

kRa 

Jm(7r  R  o/„) 

Footnote 

0.225 

0.085 

0.241 

a 

1(T4 

106 

0.225 

0.085 

0.241 

b 

0.226 

0.084 

0.241 

c 

-0.118 

-0.131 

0.176 

a 

0.1 

103 

-0.119 

-0.132 

0.178 

b 

0.136 

-0.199 

0.241 

c 

-0.192 

-0.042 

0  197 

a 

0.2 

500 

-0.193 

-0.042 

0.198 

b 

-0.213 

-0.113 

0.241 

c 

-0.025 

0.244 

0.245 

a 

0.4 

250 

-0.025 

0.245 

0.246 

b 

0.058 

0.234 

0.241 

c 

-0.413 

-0.400 

0.574 

a 

0.4975 

201 

-0.403 

-0.409 

0.574 

b 

0.241 

0.015 

0.241 

c 

“Numerical  evaluation  of  the  integral  in  Eq.  (17). 
^Evaluation  of  the  new  analytical  solution  in  Eq.  (21a). 
“Evaluation  of  Bradley  and  Hudimac’s  solution  in  Eq.  (23). 


The  function  sin  (mn9/0o)  is  plotted  for  the  first  three  modes,  corresponding  to  m  =  1,  2,  3,  in 
Fig.  2.  Note  that  each  mode  is  zero  on  the  boundaries,  and  that  the  mth  mode  exhibits  m  extrema. 

TOTAL  NUMBER  OF  MODES  IN  THE  FIELD 

The  integral  for  /„  in  Eq.  (19)  is  essentially  zero  when  v  exceeds  a  certain  value.  This  can  be 
appreciated  by  allowing  v  to  become  indefinitely  large,  in  which  case  the  trigonometric  function  oscil¬ 
lates  far  more  rapidly  than  the  rest  of  the  integrand.  If  the  slowly  varying  part,  that  is,  everything 
except  the  cosine,  is  treated  as  a  constant  and  taken  outside  the  integral,  then  the  integral  of  the  cosine 
function  is  identically  zero  (bearing  in  mind  that  v  =  mir/Oo,  and  we  are  assuming  that  00  is  a  submulti¬ 
ple  of  7 r).  The  implication  of  this  is  that  only  a  finite  number  M  of  modes  contributes  significantly  to 
the  field  at  any  point  in  the  wedge. 

To  estimate  M,  we  set 

(l (cr)  =  kR0(l  —  2a  cos  cr)1^2 

and  employ  the  condition  that  the  integral  in  Eq.  (19)  takes  appreciable  values  only  when 

mn  ,  da 
v  =  — —  ^  -F- 


3  2  1  =  m 


(24) 

(25) 
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where  (d^/do)  |max  is  the  maximum  ‘frequency’  of  the  exponential  term  in  the  integrand.  A  simple 
calculation  shows  that 


dji 

do 


max 


=  {1  -  (1  “  4fl2)12},/2- 


Thus,  the  maximum  value  of  m  from  Eq.  (25)  is 


M  =  E 


kR  o@o 


TT'fl 


{1  -  (1  -  4a2)  I/2}1/2 


(26) 


(27) 


where  the  symbol  E[  ]  denotes  "integer  part  of." 

It  is  apparent  from  Eq.  (27)  that  M,  the  total  number  of  propagating  modes  at  the  receiver, 
depends  on  the  positions  of  the  source  and  the  receiver.  There  are  two  cases  in  particular  which  are 
worth  noting. 

Case  1.  Source  and  receiver  coincident  in  r  —  z  plane  (a  =  1/2 

In  this  configuration,  we  have  r  =  r'  and  z  =  0,  giving  Rq  ==  V2  r\  and  from  Eq.  (27) 


M  =  E 


kr'90 


IT 


(28) 


kR^O^a 

—  F 

krr'OQ 

7 T 

—  L, 

7r(r2  +  r'2  +  z2)l/\ 

M  =  E 

In  certain  specific  situations,  this  may  be  simplified  further,  as  follows: 

krO  o 


Thus  the  total  number  of  modes  in  this  case  is  governed  by  the  water  depth  r'90  at  the  source. 

Case  2.  Source  and  receiver  well  separated  in  range  and/ or  z  (a  — 1 •  0) 

In  this  case,  the  right-hand  side  of  Eq.  (26)  becomes  equal  to  kR0a ,  and  the  expression  in  Eq. 
(27)  for  AT  reduces  to 

(29) 

(30a) 
(30b) 

(30c) 


and 


M  =  E 

M  =  E 

M  =  E 


7 T 

kr'6  o 


7 T 


krr'Or 


n  z 


when  r'  »  r,  z  =  0, 
when  r  »  r' ,  z  —  0, 

when  |z|  »  r  and  |z|  »  r'. 


Equations  (30a)  and  (30b)  correspond  to  the  receiver  being  upslope  and  downslope  of  the  source, 
respectively.  In  both  cases,  the  total  number  of  modes  is  governed  by  the  water  depth  (r90  or  r'60 )  at 
the  shallower  location.  Moreover,  Eqs.  (30a)  and  (30b)  are  the  same  results  that  would  have  been 
obtained  from  shallow  water  theory,  had  we  assumed  that  the  bottom  was  locally  flat  at  r  and  r',  respec¬ 
tively.  (A  similar  statement  applies  to  Eq.  (28)). 
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SPATIAL  PROPERTIES  OF  THE  MODES 

As  the  solution  for  the  acoustic  field  has  the  angular  dependence  separated  out,  the  behavior  of 
the  modes  in  the  range  and  z  directions  can  be  discussed  independently  of  their  angular  properties.  To 
examine  this  behavior,  we  should  strictly  employ  the  expression  for  /„  in  Eq.  (21a).  (In  certain  appli¬ 
cations,  for  example  in  determining  the  response  of  an  array  of  sensors  in  a  wedge-shaped  channel 
where  phase  is  of  critical  importance,  this  would  be  essential.)  However,  we  are  interested  here  only  in 
establishing  the  broad  features  of  the  field,  and  for  this  purpose  the  simpler  expression  in  Eq.  (23)  is 
adequate.  In  particular  we  shall  focus  attention  on  the  Bessel  function  Jv(kR0a),  since  this  contains 
most  of  the  interesting  properties  of  the  field;  the  exponential  functions  merely  influence  the  phase  of 
the  field,  and  R0  in  the  denominator  is  only  a  spreading  term. 

Range-Dependence  of  the  Modes  (z  =  0) 

We  begin  by  considering  upslope  propagation,  with  the  source  much  further  from  the  apex  of  the 
wedge  than  the  receiver  (i.e.,  r'  »  r ),  and  with  the  source  and  the  receiver  in  the  same  vertical  plane 
perpendicular  to  the  apex  (i.e.,  z  =  0).  Under  these  conditions,  the  Bessel  function  in  Eq.  (23)  can  be 
expressed  simply  as  Jv{kr).  This  function  is  plotted  in  Fig.  3  for  v  =  18,  36,  and  54,  corresponding  to 
the  first  three  modes  in  a  10°  wedge.  The  vertical  arrows  are  at  the  mode  cutoff  points  calculated  from 
shallow  water  theory.  The  curves  are  highly  oscillatory  beyond  some  cutoff  value  of  kr,  which  increases 
with  the  mode  number.  Below  cutoff  (i.e.,  for  smaller  values  of  kr)  the  modal  field  falls  rapidly  to  zero 
and  remains  there  right  up  to  the  apex.  Immediately  above  cutoff,  the  mode  envelope  is  a  maximum, 
and  it  decays  approximately  as  ( kr)~l/ 2  as  kr  rises. 


APEX 


TO  THE 
SOURCE 


TO  THE 
SOURCE 


TO  THE 
SOURCE 


By  examining  the  series  expansion 


= 


Or/2)1- 

vm! 


Or/2) 

(v  +  1) 


(31) 


for  the  Bessel  function  of  integer  order  v,  it  can  easily  be  seen  from  the  term  outside  the  square  brack¬ 
ets  that,  for  large  v,  Jv{x)  drops  extremely  rapidly  to  zero  when  the  argument  x  falls  below  the  order 
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v.  Taking  equality  between  x  and  v  as  the  condition  for  the  onset  of  this  behavior  (i.e.,  when 
x  v,  J„(x)  =  0),  we  can  define  for  the  m  th  mode  the  cutoff  range  from  the  apex,  rcm,  as 

krcm  =  mi t/6q.  (32) 

This  result  is  equivalent  to  that  in  Eq.  (30a),  derived  from  considerations  of  the  mode  integral. 

The  vertical  arrows  in  Fig.  3  indicate  the  cutoff  points  for  the  three  modes,  calculated  from  Eq. 
(32).  If  rcm  60  is  equated  with  the  water  depth  at  the  cutoff  range,  then  Eq.  (32)  is  the  same  as  that 
obtained  from  shallow  water  theory  for  mode  cutoff.  The  obvious  physical  interpretation  of  this  is  that 
if  the  wedge  is  regarded  as  being  locally  uniform  in  depth,  then  at  ranges  less  than  rcm  the  channel  is 
not  sufficiently  deep  to  support  the  m  th  mode. 

A  qualitative  interpretation  of  the  curves  in  Fig.  3  can  be  given  in  terms  of  rays  (Fig.  4).  A  ray 
travelling  upslope  undergoes  numerous  reflections  from  the  boundary  planes,  increasing  its  grazing 
angle  on  successive  bounces  from  either  one  of  the  planes  by  twice  the  wedge  angle.  Eventually,  as  the 
apex  is  approached,  the  grazing  angle  becomes  so  large  that  the  ray  is  turned  around  and  proceeds  to 
propagate  back  downslope.  In  the  vicinity  of  the  turnaround  point,  which  corresponds  closely  to  the 
mode  cutoff  range  given  by  Eq.  (32),  the  density  of  rays  is  high,  which  is  consistent  with  the  relatively 
high  level  of  the  modal  envelope  immediately  above  cutoff. 


Fig.  4  —  A  ray  propagating  upslope  is  eventually  turned  around 
and  then  propagates  back  downslope 


Z-Dependence  of  the  Modes 


The  argument  of  the  mode  amplitude  function  Jv(kRoa )  depends  on  z  (the  separation  of  the 
source  and  receiver  in  the  direction  parallel  to  the  apex)  through  the  term  Ro  (see  definition  in  Eq. 
(20)).  We  now  examine  this  z-dependence  downslope  of  the  source  with  the  ranges  of  the  source  and 
the  receiver  held  fixed.  Fig.  5  shows  the  function 


krr' 

(r2  +  r'2  +  z2)12 


corresponding  to  the  first  mode  in  a  10°  wedge,  plotted  as  a  function  of  z,  assuming  that  r  >>/•'.  Note 
that  when  |z|  exceeds  a  certain  critical  value  the  field  falls  to  zero,  and  where  it  is  nonzero  it  is  highly 
oscillatory  suggesting  strong  interference.  The  envelope  of  the  oscillatory  function,  derived  from  the 
asymptotic  expansion  of  the  Bessel  function,  is  shown  as  the  broken  line,  with  the  cutoff  points  in  |z| 
determined  as  in  Eq.  (32),  by  equating  the  argument  with  the  index.  The  resemblance  of  this  outline 
to  a  butterfly  with  spread  wings  is  apparent,  which  suggests  the  term  "modal  butterfly."  The  wingspan  is 
clearly  dependent  on  the  mode  number,  as  illustrated  in  Fig.  6  for  the  first  four  modal  butterflies  of  a 
10°  wedge.  As  the  mode  number  increases,  it  is  evident  that  the  width  of  the  nonzero  field  region  in 
the  z-direction  becomes  narrower,  eventually  cutting  off  altogether  above  some  value  of  m  depending 
on  the  local  conditions.  The  criterion  governing  the  number  of  propagation  modes  at  any  given  point 
has  already  been  established  and  is  given  in  Eq.  (27). 


The  significance  of  the  modal  butterflies  becomes  apparent  when  the  wingspan  is  examined  as  a 
function  of  range.  This  is  shown  for  the  first  mode  in  Fig.  7.  The  radiation  field  diverges  from  the 
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Fig.  5  —  z-dependence  of  the  first  mode  at  a  fixed  range  5  km  downsiope  of  the 
source,  which  is  located  in  the  z  =  0  plane  (0O  =10°) 


Fig.  6  —  The  first  four  modal  butterflies  at  a  fixed 
range  5  km  downsiope  of  the  source,  which  is  located 
in  the  z  —  0  plane  (0O  =  10°) 


Fig.  7  —  Range-dependence  of  the  wingspan 

downsiope  of  the  source.  Outside  the  beam  the  field 
decays  rapidly  to  form  shadow  zones  (f)0  =  10°). 
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source  to  illuminate  a  fan-like  region  beyond  which  there  is  essentially  no  energy,  and  acoustic  shadow 
zones  are  formed.  Thus,  the  interaction  of  the  radiation  with  the  boundaries  creates  a  beaming  effect. 
In  practical  situations  this  could  be  important  since  regions  of  the  ocean  are  not  being  ensonified.  Fig¬ 
ure  8  shows  the  same  phenomenon  for  the  first  four  modes,  and  clearly  illustrates  the  nested  structure 
of  the  modal  beams— as  the  mode  number  rises,  the  width  of  the  beam  decreases,  which  is  consistent 
with  the  argument  given  above  for  the  behavior  of  the  wingspan  with  mode  number. 

A  similar  situation  is  encountered  with  the  field  upslope  of  the  source,  except  that  in  this  case  the 
field  converges  as  the  receiver  recedes  from  the  source  (i.e.,  as  the  shoreline  is  approached).  This  is 
illustrated  in  Fig.  9  for  the  first  mode.  The  beam  shape  is  hyperbolic  [8],  as  is  easily  confirmed  from 
the  mode  amplitude  Bessel  function.  The  beams  corresponding  to  the  higher  order  modes  again  fall 
within  those  of  the  lower  order  modes,  although  this  is  not  shown  in  Fig.  9.  As  before,  outside  a  given 
beam  there  is  essentially  no  energy  in  the  field  associated  with  that  particular  mode.  Thus,  the  max¬ 
imum  extent  of  the  ensonified  region  in  the  \z\  direction  is  determined  by  the  lowest  order  mode  in  the 
field. 


Fig.  8  —  The  first  four  modal  beams, 
showing  their  nested  structure 
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Fig.  9  —  The  field  and  shadow  zones  of  the  first  mode  upslope  of  the  source. 
The  boundary  curve  is  a  hyperbola  (00  =  10°). 
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RAY/MODE  DUALITY 

Each  mode  in  the  wedge  is  a  superposition  of  two  components,  represented  by  the  two  terms  con¬ 
taining  the  Hankel  functions  in  Eq.  (21a).  When  the  arguments  of  these  functions  are  very  much 
greater  than  the  order  v,  the  Hankel  functions  themselves  may  be  approximated  by  their  asymptotic 
forms.  The  range  dependence  of  the  two  modal  field  terms  is  then  given  by  the  oscillatory  functions 
exp  i.-jkR\)  and  exp  (— jkR2 ),  where  from  the  definitions  in  Eqs.  (20)  and  (22),  we  see  that 

Rl  =  [(,  -  r’)2  +  z]1/2,  R2  =  [(r  +  r')2  +  z2]1/2.  (33) 

These  expressions  show  that,  in  the  ( r,z )  plane,  i?!  is  the  range  from  the  source  to  the  receiver,  and 
R  2  is  the  range  from  an  image  source  at  (—  r',0)  to  the  receiver.  Thus,  in  the  asymptotic  region  the 
two  exponential  functions  are  representative  of  two  sets  of  waves,  one  of  which  is  produced  by  the 
source  and  the  other  of  which  is  associated  with  the  image. 

The  appearance  of  the  two  field  components  in  the  modal  solution  can  be  understood  by  examin¬ 
ing  a  ray  path  in  the  wedge,  or  more  precisely,  the  horizontal  projection  of  a  ray  onto  the  sea  surface  as 
it  bounces  between  the  surface  and  the  bottom.  Weston  [9]  has  discussed  this  device  and  has  shown 
from  geometrical  arguments  that  the  shape  of  the  horizontally  projected  ray  path  is  a  hyperbola  (Fig. 
10).  That  is,  the  repeated  acoustic  interactions  with  the  inclined  planes  forming  the  wedge  introduce  a 
curvature  into  the  horizontal  direction  of  travel  of  the  ray  and  eventually  lead  to  a  turning  point  as  the 
shoreline  is  approached.  (Weston  refers  to  this  phenomenon  as  horizontal  refraction.)  If  the  asymp¬ 
tote  associated  with  the  branch  of  the  hyperbola  remote  from  the  source  is  extended  beyond  the  apex, 
it  passes  through  the  image  source  at  (—r',0). 


z 


APEX  OF 
WEDGE 


Fig.  10  —  Horizontal  projection  of  a  ray  path  in  the  wedge 


There  are  many  possible  hyperbolic  ray  paths  between  a  source  and  a  receiver  in  the  wedge.  Most 
of  these  do  not  show  constructive  interference  and  hence  do  not  correspond  to  a  mode.  The  criterion 
for  a  ray  to  correspond  to  the  m  th  mode  is  that  the  grazing  angle  of  the  ray  at  the  vertex  of  its  hyperbolic 
path  must  be  the  same  as  that  of  the  ray  corresponding  to  the  mth  mode  in  shallow  water  whose  depth  is  equal 
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to  the  depth  at  the  vertex .  Weston  [9]  has  shown  that  the  angle  between  the  asymptotes  of  a  hyperbolic 
ray  path  and  the  shoreline  is  the  same  as  the  grazing  angle  at  the  vertex.  It  follows  that  the  criterion 
given  above  for  a  ray  to  correspond  to  a  mode  fixes  the  directions  of  the  arms  of  the  modal  hyperbola. 
This  is  an  important  factor  in  interpreting  the  modal  beams  and  shadow  zones  in  the  wedge. 


The  quantity 


J* h  sin  a 
o  c 


dh, 


(34a) 


where  h  is  the  water  depth  and  a  is  the  grazing  angle,  is  a  ray  invariant  [10].  For  isovelocity  water 
(which  is  our  case),  a  and  c  are  independent  of  depth,  and  so 

T  =  -  sin  a  .  (34b) 

c 

Thus,  in  the  wedge  the  quantity  in  Eq.  (34b)  is  a  constant  along  the  hyperbolic  ray  path.  But  the  condi¬ 
tion  for  a  ray  to  correspond  to  a  mode  in  a  shallow  water  channel  with  pressure-release  boundaries  is 


h  sin  a  =  — - —  =  const.  (35) 

k 

On  comparing  these  two  expressions,  it  can  be  deduced  that  in  the  wedge  a  hyperbolic  ray  correspond¬ 
ing  to  a  mode  behaves  as  though  it  were  a  mode  in  locally  shallow  water  all  along  its  track.  That  is,  at 
each  point  along  the  hyperbolic  path,  the  upward  and  downward  traveling  wavefronts  obey  the  same 
conditions  for  constructive  interference  as  if  the  channel  were  locally  uniform  in  depth. 


For  a  fixed  source/receiver  configuration  in  the  wedge,  there  are  only  two  possible  modal  rays 
corresponding  to  the  m  th  mode.  They  are- illustrated  in  Fig.  11.  (The  parameters  of  these  hyperbolae 
are  easily  expressed  in  terms  of  r,  r',  z,  m ,  and  90  from  a  geometrical  argument,  as  detailed  in  Appen¬ 
dix  C).  One  is  a  direct  ray  from  the  source  and  the  other  is  an  indirect  ray,  which  has  been  turned 
around  on  approaching  the  apex,  as  shown  in  the  figure.  When  the  receiver  is  in  the  asymptotic 
regions  of  these  two  hyperbolic  ray  paths,  the  direct  and  indirect  rays  correspond  to  the  field  terms 
from  the  source  and  its  image,  respectively,  in  the  modal  solution  (Eq.  (21a)). 
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Fig.  11  —  Horizontal  projection  of  the  direct  and  indirect  ray  paths 
for  fixed  source  and  receiver  positions 
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We  can  now  use  the  definition  of  a  modal  ray  to  interpret  the  beaming  of  radiation  in  the  wedge, 
as  predicted  by  the  modal  solution.  Figure  12  shows  upslope  propagation  with  the  broken  line 
representing  the  edge  of  the  beam  associated  with  the  first  mode.  That  is,  the  broken  line  is  the  same 
as  the  solid  hyperbola  in  Fig.  9  obtained  from  the  modal  solution.  The  solid  lines  in  Fig.  12  represent  a 
few  of  the  first-mode  rays,  each  of  which  is  launched  from  the  source  at  a  different  azimuthal  angle. 
Notice  that  the  modal  rays  fall  precisely  within  the  beam.  Whichever  description  of  the  field  is  used, 
either  modes  or  rays,  the  energy  associated  with  a  given  mode  falls  in  the  same  well-defined  region, 
beyond  which  are  shadow  zones.  This  beaming  of  the  radiation  can  now  be  understood  from  the  ray 
argument  to  be  due  entirely  to  ray  curvature  arising  from  the  repeated  acoustic  interactions  with  the 
inclined  surface  and  bottom  planes.  The  curvature  is  such  that  no  radiation  showing  constructive 
interference  characteristic  of  the  m  th  mode  can  extend  into  the  region  beyond  the  beam  associated 
with  the  m  th  mode.  Naturally,  if  a  point  receiver  is  placed  anywhere  within  a  modal  beam,  such  as  at 
point  A  in  Fig.  12,  only  two  of  the  modal  rays  will  be  incident  upon  it,  as  indicated  by  the  double¬ 
headed  arrows.  The  beam  itself,  of  course,  is  comprised  of  many  modal  rays. 
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Fig.  12  —  Upslope  propagation  showing  the  convergent  first-mode  envelope 
and  some  first-mode  rays 


EXPERIMENTAL  EVIDENCE 

The  beaming  of  radiation  in  the  wedge  that  is  predicted  theoretically  has  been  observed  experi¬ 
mentally.  A.  B.  Wood  [11],  in  one  of  his  model  tank  experiments,  examined  the  field  in  a  wedge- 
shaped  domain.  The  tank  he  used  was  several  feet  long  with  a  depth  of  six  or  so  inches,  and  the  wedge 
was  formed  by  a  sheet  of  glass  inclined  at  an  appropriate  angle  to  the  water  surface.  To  record  the  field 
upslope  of  the  source,  he  used  a  novel  technique  which  is  appealing  in  its  simplicity.  He  first  painted 
the  glass  sheet  with  a  water  soluble  paint  (distemper),  waited  until  it  had  dried  but  not  hardened,  and 
then  placed  the  sheet  in  the  tank  to  form  the  wedge.  He  found  that  the  sound  field  left  an  impression 
in  the  paint. 

Figure  13  shows  the  result  of  one  of  A.  B.  Wood’s  experiments.  The  photograph  is  deliberately 
inverted  in  order  to  have  the  same  orientation  as  the  previous  figures,  with  the  apex  of  the  wedge  run¬ 
ning  horizontally  along  the  top  and  the  source  located  at  the  bottom.  The  interference  pattern 
displayed  by  the  field  is  easily  visible  in  the  figure,  but  the  most  prominent  feature  is  the  convergent 
beam  running  towards  the  apex.  It  is  difficult  to  make  a  quantitative  comparison  between  the  observed 
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Fig.  13  —  Acoustic  field  observed  by  A.  B.  Wood 
in  a  wedge-shaped  domain  [11]* 


beam  shape  and  that  predicted  theoretically  because  not  enough  information  is  available  on  the  experi¬ 
mental  arrangement.  The  overall  features  of  the  measured  radiation  field  are,  however,  consistent  with 
the  theory. 

CONCLUDING  REMARKS 

The  mode  and  ray  descriptions  of  the  acoustic  field  in  a  wedge  with  perfectly  reflecting  boundaries 
are  complementary.  The  modal  solution  gives  a  complete  description  of  the  field,  including  the  phase, 
and  is  thus  appropriate  to  calculations  of  spatial  coherence,  array  performance,  and  related  topics.  The 
dual  description,  in  terms  of  rays,  gives  a  pictorial  view  of  the  field  which  provides  a  valuable  intuitive 
understanding  of  propagation  phenomena  in  the  wedge.  The  formation  of  shadow  zones,  for  example, 
has  a  simple  physical  interpretation  in  terms  of  modal  rays. 

A  new  solution  for  the  field  in  the  perfect  wedge  has  been  presented  here  which  is  valid 
throughout  the  ensonified  regions.  It  is  possible  to  derive  this  solution  because  the  Helmholtz  equation 
for  the  problem  is  separable.  By  way  of  contrast,  the  Helmholtz  equation  for  the  penetrable  wedge  is 
not  separable.  This  implies  a  degree  of  complexity  beyond  that  encountered  in  the  perfect  wedge  prob¬ 
lem.  Indeed,  this  is  to  be  expected  because  when  one  boundary  is  penetrable  there  are  critical  angle 
effects  with  which  to  contend,  and  even  for  those  rays  which  undergo  total  internal  reflection,  there  is  a 
phase  change  on  reflection  from  the  boundary.  Despite  these  difficulties,  it  is  possible  to  obtain  an 
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approximate  solution  for  the  field  in  the  penetrable  wedge  by  using  the  perfect  solution  derived  here  in 
conjunction  with  an  argument  which  relies  heavily  on  the  duality  between  rays  and  modes.  This  will  be 
discussed  in  a  future  publication. 
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Appendix  A 

EVALUATION  OF  THE  NORMAL  MODE  INTEGRAL 


The  inversion  integral  giving  the  amplitudes  of  the  normal  modes  (Eq.  (19)  in  the  text)  is 

1  p*  cos  v  c t  exp  [— jkR0(l  —  2 a  cos  a-)1'2] 


/„  = 


f 

J  0 


(1  —  2 a  cos  cr)1^2 


-dcr, 


where  a  and  R0  are  as  defined  in  Eq.  (20).  Now,  with 

x  =  cos  cr, 

we  are  interested  in  the  function 

p(x)  =  (1  —  2  ax)1'2  =  (1  —  2axo)1/2 


,  a  (x0  -  x) 

1  H - - - h 

(1  -  2axn) 


(Al) 

(A2) 

(A3) 


where  the  term  on  the  right  is  the  Taylor  expansion  of  p(x)  about  x  =  xq.  On  splitting  the  range  of 
integration  into  two,  and  setting  x0  =  1  in  the  first  interval  and  x0  =  -1  in  the  second  (these  values 
correspond  to  the  turning  points  in  p(x)),  the  integral  in  Eq.  (Al)  can  be  expressed  as 


1 

ttR\ 


pir/l 

I  COS  V  <7 

Jo 


exp{— jkR\[l  +  6](1  —  cos  cr)  +  . . .  ])dcr 


where 


+ 


1  r n 

TtR2  -*n/2 


cos  v  cr 


exp{— jkR2[\  -  bi{\  +  cos  cr)  +  . . .  ]}dcr. 


(A4) 


R,  =  R0(l  -  2a)1'2,  R2  -  R0(l  +  2a)1/2  (A5a) 

bx  =  a/{\  -  la),  Z>2=  a/d  +  2a)  (A5b) 


and  p(x)  in  the  denominator  has  been  approximated  by  the  first  term  in  the  expansion  in  Eq.  (A3). 
The  two  components  of  Iv  can  now  be  written  in  the  form 


where 


and 


exp  K/fcftjXl  +  fri)]  exp  [-jkR2(  1  -  b2)\ 
7 tR\  1  7TR2 

/•  ir/l 

I\  =  cos  v  c t  exp  [jikRibi  cos  o-)]do- 


J*  7 T 

cos  v  a  exp  [y(fcR2*2  cos  c r))d<T. 

7T/2 


(A6) 

(A7a) 

(A7b) 


In  these  expressions,  the  terms  beyond  those  shown  in  the  expansions  in  Eq.  (A4)  have  been 
neglected.  This  is  justified  because  the  major  contributions  to  the  integrals  in  Eq.  (A4)  come  from 
around  the  turning  points  of  p(x),  where  only  those  terms  shown  are  significant. 


With  a  little  algebraic  manipulation,  the  integrals  in  Eq.  (A7)  can  be  expressed  as  standard  forms 
which  can  be  found  in  most  tables  of  integrals  (e.g.,  Ref.  12).  For  v  an  integer,  the  results  are 


and 


h 


{kRxbx) 


(A8a) 
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H™  ( kR2b2 ), 


(A8b) 
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which,  when  substituted  into  Eq.  (A6)  give 


T  1  .VI T 


exp  1-jkRi (1  +  bx)]H^_  (kRxbx)  exp  [-yfcfi2(l  ~  b2)]H™  ( kR2b2 ) 


Ri 


R 7 


This  is  the  result  in  Eq.  (21a)  in  the  text  that  we  set  out  to  prove. 


(A9) 
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Appendix  B 

THE  FIELD  IN  THE  VICINITY  OF  THE  SOURCE 

The  field  from  a  point  source  within  a  bounded  domain  shows  a  singularity  at  the  source  point, 
and  in  the  immediate  vicinity  of  this  point  it  has  the  same  form  as  the  free  field  of  the  source.  Indeed, 
this  fact  is  often  employed  to  construct  solutions  for  boundary  value  problems:  the  field  is  represented 
as  a  superposition  of  the  free-field  component  and  another  component  chosen  to  satisfy  the  boundary 
conditions. 


In  the  case  of  the  wedge,  the  solution  for  the  field  in  Eqs.  (15)  and  (21a)  has  been  derived  by 
using  transform  methods,  rather  than  the  superposition  technique,  but  it  should  nevertheless  be  valid 
in  the  immediate  vicinity  of  the  source  point.  That  is,  the  modal  sum  in  this  region  should  reduce  to 
the  free-field  form.  We  show  here  that  this  is  the  case. 


When  the  receiver  is  very  close  to  the  source,  only  the  term  containing  the  Hankel  function  of 
the  first  kind  in  Eq.  (21a)  is  significant.  This  Hankel  function  has  a  large  argument  (kRyby)  and  a  large 
order  v.  Bearing  in  mind  that  at  the  source  point  there  is  a  finite  number  of  propagating  modes  (see 
Eq.  (28)  in  the  text),  it  is  always  possible  to  choose  kRyby  >  v.  (kRyby  increases  as  the  source  point 
is  approached).  On  setting 


kRyby  =  v  sec  p, 

the  Hankel  function  of  the  first  kind  can  be  expressed  approximately  (see  Ref.  Bl) 

(kRyby)  =  H™  (v  sec  p)  -  _  /—  2 

V  7r  v  tan  p 

From  the  definition  of  p  in  Eq.  (Bl),  we  have 


as 


exp 


^(tan  b  —  p)  — 


IcRyby 

tan  p  —  p  =  — 7t/2  H - 1- 


IkRyby 


+  ...  . 


Assuming  that  v/(kRyby)  «  1  for  all  v,  so  that  the  series  may  be  truncated  as  shown,  this  gives 


M,b0  “  V'^tibrexp ^ _  “,'’1 


— - - 1-  — 

IkRyby  4 


Thus  the  expression  for  Iv  in  Eq.  (21a)  can  be  approximated 


as 


exp  (—j  7t/4)  exp  '  v 

Al 


R\ 

When  this  substituted  into  Eq.  (15),  the  expression  for  the  field  becomes 

J 

e 


IkRyby 


<f>  =  ^  (27 rkRyby)  1/2  — —  ^ — —  exp  (-j  n/4 )  ^  sin  v  6  sin  v  O'  exp 

1  u 


2kR,b, 


(Bl) 


(B2) 


(B3) 


(B4) 


(B5) 


(B6) 


To  evaluate  the  summation  in  Eq.  (B6),  we  first  employ  the  trigonometric  identity 

sin  v  9  sin  v  0'  =  1/2 [cos  v(9  -  O')  -  cos  v(0  +  9')],  (B7) 
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but  retain  only  the  difference  term,  since  the  remaining  term  makes  a  negligible  contribution  to  the 
field  when  9  — *■  9'.  The  sum  in  Eq.  (B6)  can  therefore  be  written  in  the  form 


1 


S  =  —  £  cos  v(0  —  9')  exp 


2kR\b\ 


By  expressing  the  cosine  function  as  the  sum  of  two  exponentials,  this  sum  can  be  written  as 


S  =  -J-  (5+  +  Sj. 
4 


where 


M 


S±  -  £  exp 

m— 1 


J 


±  ~~~  (9  —  9')  +  mV 


#0 


29q  kR  ]  b\ 


(B8) 

(B9) 

(BIO) 


Here  we  have  substituted  v  =  mir/90,  and  the  upper  limit  on  the  summation  is  given  by  Eq.  (28). 
Now,  provided  the  maximum  value  of  the  term  containing  m 2  in  the  argument  of  the  exponential  func¬ 
tion  in  Eq.  (BIO)  is  greater  than  unity,  S±  may  be  approximated  by  the  following  integral: 


S+  =  S_  =  -L  f 

+  Vu  J  o 


exp 


>^2  + 


2  wy 


dy. 


where 


u  = 


7 T 


.  ■,  2t»-f  (8-  9'). 

efakR^by  90 


(Bll) 


(B12) 


The  integral  in  Eq.  (Bll)  is  a  standard  form  [B2]  which  can  be  expressed  exactly  in  terms  of  Fresnel 
integrals.  However,  in  the  limit  as  9  —'■9'  (which  is  our  case)  the  contribution  from  the  Fresnel 
integrals  falls  to  zero,  and  we  find  that 


2  S 


e  „  1  fir  .  w2 

=  S+  =  S_=  — _ /  —  exp  -y — 
2  V  u  u 


exp  O'  if! 4). 


(B13) 


When  this  expression  for  S  is  substituted  into  Eq.  (B6),  in  place  of  the  summation,  the  field  in  the 
vicinity  of  the  source  is  found  to  be 


Q  exp  (-jkR{) 
- R, - exp 


-J- 


u 


The  argument  of  the  second  exponential  function  in  this  expression  is  now 

.,2 


t  (9-  9')1 


u 


=  kRxb 


=  kR\b\ 


1  - 

.  (9  -  9')2 1 

k 

2 

(B14) 


(B15) 


a*  kRxbi  11  -  cos (9  -  9')], 

where  the  approximation  holds  because  we  are  specifically  interested  in  the  case  where  9  — 1 ►  9'.  Now,  if 
R  is  the  distance  between  the  source  and  the  receiver,  then 

J?2  =  r2  +  r'2  +  z2  —  2 rr'  cos  (9  —  9')  (B16) 

=  R 2  +  2rr'[\  -  cos  (9  -  9')] 

so  that 

>2_p.2  (R  _  Rl)(R  +  Rl) 


[1  -  COS  (9  -  flOl  = 

i\  ^  Zi\  i 


22?, 


(B17) 


=  (R  -  Rx) 
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where  the  approximation  is  valid  because  in  the  immediate  vicinity  of  the  source  R  =  R{.  Since 
kRxb\  =  krr'/Ri,  it  follows  from  Eqs.  (B15)  and  (B17)  that 


(B18) 


and  hence  the  expression  for  the  field  in  Eq.  (B14)  takes  the  form 

.  =  _Q_  exp  (- jkR ) 

9  4tt  R 


(B19) 


where  the  denominator  has  been  written  as  R  rather  than  Rv  But  the  expression  in  Eq.  (B19)  is  pre¬ 
cisely  the  free-field  produced  by  the  source,  which  is  what  we  set  out  to  prove.  Note  that  Bradley  and 
Hudimac’s  [B3]  solution  for  the  field  in  the  wedge  (Eq.  (23))  does  not  reduce  to  the  required  form  in 
the  vicinity  of  the  source  point. 
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Appendix  C 

MODE/RAY  HYPERBOLAE 


The  two  hyperbolic  ray  paths  associated  with  a  particular  mode  are  illustrated  in  Fig.  Cl.  The 
parameters  of  the  two  curves  are  derived  below. 
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Fig.  Cl  —  Horizontal  projection  of  the  direct  and  indirect  ray  paths 
for  fixed  source  and  receiver  positions 


Let  the  point  (r,z)  fall  on  one  of  the  hyperbolae.  The  Eq.  for  the  hyperbola  can  then  be  written 
as 


f_  _  (f  -  *m)2  =  . 

d  bl 

where  rm,  bm,  and  zm  are  the  parameters  of  the  hyperbolae  associated  with  the  m  th  mode, 
parameters  take  dual  values,  corresponding  to  the  two  ray  paths. 


(Cl) 

These 


From  shallow  water  theory,  the  grazing  angle  am  of  the  m  th  mode  at  the  turning  point  is  given 
by 


sin  am 


rtnr 

krm00 


(C2) 


where  rmQ0  is  just  the  depth  at  the  vertex  of  the  hyperbola,  and  qm  =  m7r/(k60).  Now,  the  angle 
between  the  apex  of  the  wedge  and  the  asymptotes  of  the  hyperbola  is  also  equal  to  am  [Cl],  and  hence 
we  have 
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It  follows  from  Eqs.  (C2)  and  (C3)  that 


tan  «m  =  -f-. 

K 


(C3) 


^)1/2- 


(C4) 


We  now  require  two  more  conditions  to  determine  the  three  parameters  of  the  hyperbola  They  are 
that  the  curve  must  pass  through  the  source  point  (/•', 0)  and  the  receiver  point  ( r,z )  (both  rotated  onto 
the  surface).  These  conditions  lead  to  the  equations 


and 


(C5) 


r2 


(z  ~  zj2 


=  1. 


rm  bl 

On  eliminating  zm  from  these  equations,  we  find  that 

2  -  -  O1'2  ±  <r°  -  r2)'«](r2  -  » 

Assuming  that  r2  »  »  ri  and  ,  su  then  Eq.  (c7)  may  be  eas„y  solved  to  giye 


rm  Qn 


1  + 


1/2 


( r  ±  r’)2 

Simple  algebra  shows  that  the  remaining  two  parameters  are  given  by  the  expressions 

*2 


zm  =  ± 


ra~  q1 

‘  H  m 


l  + 


(r  ±  r')2 


(r  +  r ') 


and 


(C6) 

(C7) 

(C8) 


(C9) 


bm  = 


1  + 


(r  ±  r')2 


(r  ±  r') 


(CIO) 


Where  there  is  a  choice  of  sign  in  Eqs.  (C8)  to  (CIO),  which  are  the  results  we 
relates  to  the  indirect  ray  path  and  the  minus  sign  to  the  direct  ray  path. 


require,  the  plus  sign 
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